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Abstract 

We study low-lying states of even-even carbon isotopes in the range A = 10 — 20 within the 
3 ■ large-scale ah initio no-core shell model (NCSM). Using several accurate nucleon-nucleon (NN) as 

well as AW plus three-nucleon (NNN) interactions, we calculate excitation energies of the lowest 
J> ' 2^ state, the electromagnetic B(E2; 2^ — t- 0J~) transition rates, the 2^ quadrupole moments as 

, measure 2"'"-state lifetimes indicate an interesting evolution of nuclear structure that is a challenge 

to reproduce theoretically from first principles. Our calculations do not include any effective charges 
or other fitting parameters. Overall, we find a good agreement with the experimentally observed 
. trends, although our calculated B(E2; 2^ — > 0^) value for ^^C is lower compared to the most recent 

, measurements. Relative transition strengths from higher excited states are investigated and the 

influence of NNN forces is discussed. In particular for ^^C we find a remarkable sensitivity of the 
transition rates from higher excited states to the details of the nuclear interactions. 
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I. INTRODUCTION 



Electric quadrupole (E2) matrix elements are important quantities in probing nuclear 
structure. In particular, they are very sensitive to nuclear deformation, the decoupling of 
proton and neutron degrees of freedom, and they are often affected by small components of 
the nuclear wave function. In this paper we study systematically the observables that are 
obtained from diagonal and non-diagonal E2 matrix elements in the range of even carbon 
isotopes, from to the very neutron rich ^'^C. Quadrupole moments, corresponding to 
diagonal E2 matrix elements, are inherently difficult to measure for excited 2+ states. Off- 
diagonal matrix elements, however, have recently been studied for several heavy carbon 
isotopes using lifetime measurements [l-5|. In this way, the reduced transition probability, 
B (E2; 2^ — )■ 0^), can be extracted since it's inversely proportional to the lifetime of the 
2^ state. As a result of these experimental studies several different claims have been made 
on the nuclear structure in this chain of isotopes. Initial excitement was triggered by the 
observation of a strongly quenched E2 transition in ^^C Based on the liquid-drop model, 
which predicts the i?(E2) to be inversely proportional to the 2"*" excitation energy, Imai et 
al. jl| claimed an anomalous reduction of the E2 strength when comparing 2"*" lifetimes for 
(^2+ = 7.01 MeVland ^^C (^2+ = 1-77 MeV). The ^^C{2+) lifetime was remeasured 
by Wiedeking et al. j2| providing a much shorter value, thus indicating a larger i?(E2) 
strength. Their results were analyzed in terms of shell model calculations. Adjusting the 
effective neutron charge to reproduce their measured lifetimes they made the claim that the 
results for ^^C are "normal" to this region. Lifetime measurements of i6,i8q ^g^g reported 
by Ong et al. j^. The presented results for ^^C came from a reanalysis of the original 
data ij, now giving a larger but still quenched B{E2) strength, while the new data 
indicated the persistence of the quenching of E2 strengths in heavy carbon isotopes. Possible 
explanations were put forward in terms of the decoupling of protons and neutrons resulting 
in very low values for the neutron effective charges and/or the appearance of a new proton 
magic number Z = 6 in this region. Some of these statements were backed up by new shell 
model calculations by Fujii et al. reproducing the ^^.isq rgg^its employing exceptionally 
small effective charges. An explanation in terms of core polarization effects was recently 



proposed by Ma et al. 



7|. They used a microscopic particle- vibration approach to compute 



core polarization effects on valence nucleons. In contrast with empirical effective charges. 
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usually employed in shell-model calculations, they noted a very strong quenching from core 
polarization on sd-shell neutrons for heavy carbon isotopes. 

These developments provide a strong motivation to perform ab initio calculations, without 
fitting parameters, to study the evolving nuclear structure in the carbon chain of (even) 
isotopes with particular focus on 2^ states and quadrupole moments. We have, therefore, 
carried out large-scale ab initio no-core shell model (NCSM) calculations for low-lying 
states of the even-even carbon isotopes with A = 10 — 20. These calculations are performed 
starting from realistic Hamiltonians without adjustable parameters. In particular, since our 
many-body scheme does not involve an inert-core approximation we will use bare charges 
when evaluating electromagnetic observables. Particular efforts will be made to quantify the 
uncertainties of the observables. 



A. Theoretical formalism 

The NCSM method has been described in great detail in several papers, see e.g., Ref. joj. 

Here, we just outline the approach as it is applied in the present study. We start from the 

intrinsic Hamiltonian for the A-nucleon system Ha = Trd + V, where T^ei is the relative 

kinetic energy and V is the sum of nuclear and Coulomb interactions. The potential term 

will always contain two-body operators, but we can also include three-body terms originating 

from an initial NNN force or induced by a unitary transformation of the Hamiltonian. This 

transformation is employed to soften the Hamiltonian for the description of many-body, 

low-energy observables and will be described below. 

In this work we have used several different nuclear Hamiltonians: the pure AW inter- 
1 — ' n — 

actions CD-Bonn 2000 [U^ (CDB2k), based on one-boson exchange theory, and INOY 11 1 
that introduces a nonlocality to include some effects of three-nucleon forces and is fitted 
also to three-nucleon obervables. In addition we have used the most recent chiral NN plus 
NNN interaction, i.e. the N^LO NN interaction of Ref. (jj] and a local chiral N^LO NNN 



potential with low-energy constants determined entirely in the three-nucleon system 13|. 
All three Hamiltonians reproduce NN phase shifts with very high precision. 

We solve the many-body problem in a large but finite harmonic-oscillator (HO) basis 
truncated by a maximal total HO energy of the A-nucleon system. The many-body model 
space is usually characterized by the cutoff-parameter N^^^, giving the maximum number of 
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HO excitations above the unperturbed A-nucleon ground state. The diagonahzation of the 
Hamiltonian in this many-body basis is a highly non-trivial problem because of the very large 
dimensions that we encounter. To solve this problem, we have used a specialized version 
of the shell model code Antoine jl^, adapted to the NCSM [isl. For the runs involving 
explicit NNN interactions we used the NSuite package 16l,[l7], which is in addition capable 
of performing the importance-truncated NCSM calculations described below, and also the 
code NCSD 



Due to the strong short-range correlations generated by the AW potentials, we must 
calculate an effective interaction to speed up the convergence. We employ two different 
similarity transformations to construct the effective interactions: For CDB2k and INOY 
as initial NN interactions we derive two-body effective interactions appropriate to the low- 
energy basis truncation by a unitary transformation in the two-nucleon HO basis (Lee- 
Suzuki effective interaction js], [lol). For the chiral NN + NNN Hamiltonian we employ 
the similarity renormalization group (SRG) with the initial and induced three-body terms 
included consistently 17, Q. 

Our results exhibit dependence on Wmax and the HO frequency Q that should disappear 
once a complete convergence is reached. This implies that W^ax- sequences obtained at 
different HO frequencies should all converge to the same result. This feature can be utilized 



to perform a constrained fit to multiple sequences (2l|, |22|. To this end, we use as large 
an Wmax basis as feasible for a wide range of HO frequencies, and extrapolate calculated 
observables to infinite space. Results obtained for a range of frequencies were used in the fits. 
We found that the convergence with increasing Wmax for energy observables can be well fitted 
by: X = Xoo+Co exp(-CiiVmax)/A^max, and for EM observables by x = Xoo+Co/Wmax+Ci/Wmax- 
The parameters cq and ci are allowed to vary for each Wmax- sequence, while the parameter 
Xoo is common to all Wmax- sequences and gives the extrapolated result at W^ax oo. An 
error estimate is made based on repeating the constrained fit keeping either the three highest 
or the three lowest frequencies in the selected range. 

The current limit on the fuU-Wmax-space calculations with a two-body Hamiltonian is, 
e.g., Wjnax = 6 for ^^C with a dimension of 1.4 x 10^ . To reach still larger Wmax we employ 
the importance-truncated (IT) NCSM scheme 16 



23l |. It makes use of the fact that many of 



the basis states are irrelevant for the description of a set of low-lying states. Based on many- 
body perturbation theory, one can define a measure for the importance of individual basis 
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states and discard states with an importance measure below a threshold value, thus reducing 
the dimension of the matrix eigenvalue problem. Through a sequence of IT calculations 
for different thresholds and an a posteriori extrapolation of all observables to vanishing 
threshold, we can recover the full NCSM results up to extrapolation errors jl6|. For i^.isq 
the IT scheme allows us to extend our calculations to A^max = 8 and to improve the reliability 
of our extrapolations. 



II. RESULTS 



A. Convergence and error estimates 

For our detailed studies of observables we are looking for the region in which the A^max- 
convergence is the fastest and the dependence on HQ is the smallest. This optimal frequency 
range can vary between different observables and different isotopes. We will use the A^max- 
dependence of the binding energy and the first excitation energy in the largest model 
spaces as our primary criterion for selecting the optimal frequency range. The largest model 
spaces that we were able to reach in the full-space NCSM calculations spanned from A^max = 
10 in lOQ, via A^^^x = 8 in ^^-i^c, to N^^^ = 6 in ^^'^^C and N^^^ = 4 in ^oc. The largest 
matrix dimension was D = 1.4 x 10^ for ^®C. However, using the IT-NCSM scheme we are 
able to obtain results also with A'^max = 8 for ^^'^*C. From plots such as Fig. [T|, focusing in 
particular on the trend for large model spaces, we find that the hQ-Ta.nge 10-14 MeV seems 
to work very well for these observables. This choice applies to the CDB2k-interaction for 
the whole range of carbon isotope. 

Figure [2] shows several examples of the constrained-fit procedure for quadrupole moments 
and B{E2) strengths. The results are plotted as a function of 1/A^max for the selected 
range of HO frequencies. Infinite model space corresponds to 1/A^max ^ 0. The grey error 
bands indicate the bounds obtained when fitting to A'^max-sequences for different hQ ranges. 
They form the basis for the uncertainty estimates. We note that the use of a range of 
frequencies usually include sequences that converge from above and from below. This allows 
a more precise determination of the extrapolated, final result. A particular exception to this 
behavior is the B{E2) strength of ^^C, for which all sequences converge from below. This 
will be further commented below. 
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FIG. 1: (Color online) /tO-dependence for the ground-state energy (presented as Egg/A) and the 
first excitation energy for ^^'^^C. Each curve corresponds to a particular model space represented 
by the parameter A^max (see text for details). 
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FIG. 2: (Color online) Model-space dependence of calculated E2 observables for 12,16,18q |^]^g 
NCSM. Results obtained with the CDB2k AW potential are presented as a function of l/Wmax- 
Filled (open) symbols correspond to full (importance-truncated) space results. Dotted lines corre- 
spond to constrained fits to Wmax-sequences at a fixed HO frequency. See also Table HI 

B. Binding energies and 2+ state properties 

In Figs. [3] and m we compare the calculated and experimental trends for binding energies, 
2f excitation energies, and E2 observables for the carbon chain of isotopes. It is clear that 
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TABLE I: NCSM calculated E2 observables of ^"~^^C compared with experimental results. NCSM 
results are obtained using the CDB2k interaction. The recommended values come from constrained 
extrapolations (see text for details). 





£;exc (2^) [MeV] 


Q {2t) [efm'] 


B{E2; 2+ ^0+)[e2fm^] 


12c 
14c 
16c 
18c 
20c 


3.62(2) 3.354 
3.49(3) 4.439 
5.36(72) 7.012(4) 
2.25(23) 1.766(10) 
1.83(9) 1.620(20) 
1.72(7) 1.588(22) 


-1.1 ±1.2'' — 
+6.19(20) +6(3) 
+4.71(23) — 
-3.21(19) — 
+3.82(4) — 
+4.20(32) — 


10 ±2" 12.4(2.0) 
8.83(89) 7.94(66) 
5.30(46) 3.74(50) 
2.15(9) 2.6(9)^ 4.15(73)'^ 
4.22(32) 4.3(1. 2)^^ 
4.68(66) < 5.7^ J.btli" 



"Strong mixing of tlie first two 2+ states. The estimate of E2 observables is obtained by studying the 

sums and ratjols of results for both 2+ states. 
''From Ref. 
''From Ref. 
^From Ref. U 
''From Ref. [5| 

the CDB2k interaction underbinds these isotopes by 10-20 % while the INOY interaction 
provides additional binding. The positive two-neutron separation energy for ^^C is not 
reproduced, which is probably associated with a slower convergence of states with neutrons 
in the sd shell. Excitation energies are much easier to converge and we find a very good 
agreement with the experimental trend. The possible exception is the large 2^ excitation 
energy of ^^C that is overpredicted with the INOY interaction. 

Numerical results for 2f excitation energies and E2 observables are presented in Table [H 
We stress that no effective charges were used in our NCSM calculations and that there are 
no adjustable parameters. Our calculated B (E2; 2f — > Of) agrees rather well with the most 
recent experimental data for the entire chain of isotopes. We note that the extrapolation 
of our ^^C B (E2; 2f Of) results is challenging as, unlike for other carbon isotopes, the 
B{E2) value increases with A'^max in the whole investigated HO frequency range, see Fig|2j 
This makes the extrapolation procedure less accurate. Although our final result agrees well 
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FIG. 3: (Color online) NCSM calculated binding energies and 2"*" excitation energies of 
compared with experimental results. Theoretical NCSM results are constrained extrapolations 
(see text for details). 
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FIG. 4: (Color online) NCSM calculated E2 observables of io~20q compared with experimental 
results. See also Table [H NCSM results are obtained from constrained extrapolations (see text for 
details) . 
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with the data from RIKEN [3| it is shghtly below the most recent experimental result from 
LBNL {2!. In addition, we note that our calculated quadrupole moment for the first 2"*" state 
of ^^C is Q = —3.21 ± 0.19 efm^ while for A = 12, 14, 18, 20 we find the quadrupole moment 
of the 2f state to be positive. 

A qualitative understanding of these findings can be obtained by studying the mean 
occupation numbers of different single-particle states in the NCSM wave functions. In Fig. [5] 
these occupancies are plotted for the ground- and first 2"''-state in the whole range of carbon 
isotopes. In particular for I'^'I^q ^j^g excitation mechanisms are quite obvious. In ^^C the 
2^ state corresponds to a proton excitation within the p shell, while in ^^C the 2^ state is 
obtained through a re-configuration of neutrons in the sd shell. The value of the B{E2) is 
quite sensitive to the fine details of the re- configuration and, as it involves the sd shell, we 
have a slower convergence in the NCSM. 

For we observe a very strong mixing of the first two 2"*" states using the CDB2k 
interaction at small frequencies. To get at least crude estimates of the E2 properties of the 
2f state we used a slightly different extrapolation approach: The ratios of, e.g., Q{2f) and 
(5(2^) was plotted for larger frequencies where the mixing is not observed, while the sum 
was plotted for the full range of frequencies. From such plots, for Q and i?(E2) observables, 
we can deduce estimates for Q{2f) and B (E2; 2^ — > 0^) and their uncertainties. These are 
included in Table |T] and Fig. HI 

We note that the different AW interactions used in this study give very similar isotopic 
trends for E2 observables, but with a consistently smaller magnitude for the INOY interac- 



ous. 



24 



25|. 



arge nuclear density generated by 



tion. This observation is connected to the anoma 
this interaction found already in ^He calculations 

Finally, a study of the characteristics of the second 2+ state in these isotopes strengthens 
the conclusion of the prominence of ^^C in the structural evolution of the chain of even 
carbon isotopes. The sign of the quadrupole moment of this state, (^(2^), is reversed from 
Q{2^). I.e., it's negative for all isotopes except for ^^C. In addition, as summarized in 
Table [Tll the relative B{E2) strength from this state to the ground state is much smaller 
than from the first 2^ for all isotopes but ^^C. These findings are obtained with both A7V 
Hamiltonians used in this study. However, the relative transitions from the second 2"^ in 
i6,i8q gtand out with clear differences in the predictions of CDB2k and INOY, see Table HTl 
Note, however, that the convergence of the second 2+ state is computationally more difficult. 
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FIG. 5: (Color online) Occupation numbers for the ground- and first 2"'"-state in io-20q obtained 
with the CDB2k interaction. Black (red) bars on the left (right) correspond to proton (neutron) 
occupation numbers. Filled (open) bars correspond to the ground (2+) state. 

and therefore the statements on relative transition strengths are based on runs performed 
at a single HO frequency. For ^®C, in particular, there is a strong /if2-dependence for the 
INOY results that makes the corresponding claim of a strong 2^ 2f E2 transition less 
robust. For ^^C, however, the interaction dependence is solid and intriguing. As the INOY 
interaction often hints to possible structural influence from NNN forces we continue our 



12 



TABLE II: Relative B(E2) values for transitions among excited states of ^''C. Results obtained 
at fixed HO-frequency with the CDB2k {hn = 12 MeV) and INOY {HQ = 17-18 MeV) NN 
interactions are compared. 





B{E2;2+ ^0+) 
B{E2; 2+->0+) 






B{E2;2+ ^2+) 
B(E2; 2+^0+) 




A 


CDB2k 


INOY 


CDB2k 




INOY 


14 


0.001 


0.000 


0.48 




0.81 


16 


2.2 


0.30 


2.0 




0.79 


18 


0.018 


0.22"^ 


0.047 




1.7'^ 


20 


0.017 


0.035 


0.12 




0.28 



Tor this particular interaction we observe considerable mixing between two 2+ states, with different 
structure, for certain choices of the HO frequency. These results are for hfl — 18 MeV. 



study in the next section with a more detailed investigation of the ^^C structure using chiral 
NN+NNN Hamiltonians. 



C. Higher-lying states of ^^C and the role of the NNN interaction 



Transitions from higher excited states of ^^C were also studied in a recent experiment [26|. 
In particular, the transitions 2^ ^ 2^^, 4^'" — t- 2^ and 3^^ 2^ were observed. Interestingly, 
no transition from the 2^ state to the ground state was seen. We performed additional 
calculations with different Hamiltonians to study higher excited states in ^^C and their 
electromagnetic transitions. In Fig. [6l we show the calculated and experimental energy 
levels of ^^C, and in Table [XTTl we summarize our calculated B{E2) values among excited 
states normalized to B (E2; 2f —t- Of). In particular, we compare results obtained with 
SRG-transformed chiral AW and chiral NN+NNN interactions, including the SRG-induced 
three- nucleon terms in both cases as discussed in Ref . , to those obtained with the CDB2k 
interaction. A striking feature is a strong suppression of the 2^ 0^ transition when the 
initial NNN interaction is included. Clearly, the calculation without the NNN interaction 
contradicts the new MSU experiment j26|. From Table IIIII we observe that relative E2 
transition strengths obtained with the chiral NN interaction are similar to the ones obtained 
with the CDB2k interaction. Furthermore, we see from Table HIl that the relative i?(E2) 

13 



% 4 



o 

• ^ 

• ^ 

(>< 



3- 







CDB2k 


chiral NN 


chiral NN+NNN 












' ' 2+- 






— 0+- 








% 












II 


II 


■ 1 ■ 1 ■ 1 0+- 



8 Exp. 2 4 6 Exp. 2 4 6 Exp. 



max 



max 



max 



FIG. 6: (Color online) Excitation energies of the lowest states of ^^C. Calculations using the Lee- 
Suzuki-transformed CDB2k potential at hQ = 12 MeV (left) and the SRG-evolved chiral AW and 
NN+NNN interactions with A = 1.88 fm~^ for hO, = 16 MeV (middle and right) are compared for 
different values of Wmax to experiment. 

calculated with the INOY interaction (that mimics some NNN effects) resemble results of 
the chiral NN+NNN Hamiltonian. The sensitivity to the presence of the NNN interaction is 
remarkable. The 2^ -+ Of transition is suppressed by a factor of ~ 20 in the calculation with 
the NNN. The excitation energies of the lowest five ^^C excited states are also influenced 
by the NNN interaction as seen in Fig. |6l The agreement with experimental spectrum is 
quite reasonable in all presented cases, although slightly improved in the calculation with 
the chiral NN+NNN Hamiltonian. 

From Table UTTl we also note a strong sensitivity of the 3^ -+ 2f transition to the presence 
of the NNN interaction. The calculation with the chiral NN+NNN Hamiltonian predicts a 
strongly suppressed i?(E2; 3^^ -+ 2f) transition. This transition is observed, however 26 |. 
Our calculation with the NNN interaction predicts this transition to be of Ml character as 
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TABLE III: Relative S(E2) values for transitions among excited states of ^^C. Results obtained 
with the CDB2k AW potential, the chiral NN, and the chiral NN+NNN interaction are compared. 
For CDB2k the Lee-Suzuki effective interactions were used {hO, = 12 MeV, Amax=6) and for the 
chiral interactions the SRG-evolved interactions (A = 1.88 fm~^, HQ = 16 MeV, Amax=6) including 
the induced three-nucleon terms were used. 

CDB2k chiral NN chiral NN+NNN 

2l ^Of 1 1 1 

of 2.2 0.75 0.11 

2+ 2+ 2.0 1.7 0.65 

3^^ ^ 2+ 0.36 0.31 0.02 

4+ ^ 2+ 0.89 0.69 0.80 



TABLE IV: Absolute values for magnetic dipole moments and B(M1) transition strengths of ex- 
cited states in ^•^•^^C. Results obtained with the CDB2k NN potential and the chiral NN+NNN 
interaction are compared. B{M1) in fij^ and fj. in fi]\f. Parameters as in Table Hill with Amax=4 for 



C. The brackets indicate the uncertainties of the threshold extrapolation for the IT-NCSM. 







16c 




20c 






chiral 




chiral 




CDB2k 


NN+NNN 


CDB2k 


NN+NNN 


B(M1; 2+ ^ 2+) 


0.013 


0.063 


0.015 




B(M1; 2+) 


0.17 


0.17 


0.013 






0.13 


-0.42 


0.22 


0.001(8) 


^(2^) 


1.3 


-0.79 


0.58 




^(3^) 


-3.2 


-3.1 


0.016 





seen from Table llVi We also observe a sign change of the magnetic moments of both the 2^ 
and the 2^ states in calculations with the NNN interaction included. The magnetic moment 
of the 3i state is unaffected, however. The sensitivity of the 2^ magnetic moment to the 
NNN interaction we also find in (see Table \W\ . 

Overall, we find a strong sensitivity of the electromagnetic observables in ^^C to the 
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details of nuclear Hamiltonian. More detailed experimental study of higher excited states 
and their transitions would be quite useful. 



III. CONCLUSION 

In summary we have computed low-lying states of even carbon isotopes with A = 10 — 20 
within the ab initio NCSM. We used several accurate nucleon-nucleon (AW) as well as NN 
plus NNN interactions and calculated excitation energies of the lowest 2+ state, the electro- 
magnetic B (E2; 2^ — 7- 0^) transition strengths, the 2f quadrupole moments and selected 
electromagnetic transitions among higher excited states. Unlike in the phenomenological 
shell model, our calculations do not include effective charges or any other fitting parame- 
ters. 

Overall, we find a consistent NCSM description of the B (E2; 2f — 0^^) dependence on 
the mass number for the whole carbon isotopic chain from A = 10 to 20. However, our 
calculated B (E2; 2f Of^ values for ^^C, with different Hamiltonians, all underestimate 
the most recent experimental measurements. Our results are quite similar to the findings 



of Ma et al. j?], who used a microscopic particle- vibration approach to compute core polar- 
ization effects. In their picture the reduced i?(E2) strength in heavy carbon isotopes can 
be traced back in particular to a strong quenching from core polarization on sd-shell neu- 
trons. In our approach, however, there is no such separation into core and valence degrees 
of freedom. 

In addition, we find a remarkable sensitivity of the transition rates from higher excited 
states in ^^C to the details of the nuclear interactions. The chiral NN + NNN interaction 
gives the excitation spectrum of ^^C in a slightly better agreement with experiment than the 
CDB2k NN potential and, further, the former interaction predicts the suppression of the 
2^— >0]^ transition in agreement with experimental observations. We find a strong sensitivity 
of the magnetic moments of the 2f state to the nuclear interaction in ^^C and ^''C and even 
more so for the 2^ state in ^^C. 

The NCSM calculations predict sign changes of the 2f quadrupole moments between 
different carbon isotopes. In particular, we predict a negative quadrupole moment in ^^C, a 
very small quadrupole moment in ^°C and a B (E2; 2f ^ Oj*") value in ^°C that is about the 
same as that in ^°Be. In ^^C, we obtain Q(2f)=+6.2{2) efm^. It will be worth measuring 
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these moments in the future. 

Finally, we note that the truncation of the many-body basis used in the NCSM should 
in principle be followed by a transformation of the transition operator that is consistent 
with the renormalization of the Hamiltonian. For long-range operators, such as E2, this 
transformation is not expected to produce a very different end result for calculated observ- 
ables 
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I . In addition, the small uncertainty associated with the approximation of using 
bare operators is partly built into the error estimates that we obtain from using several 
values of hfl and A/^max- 
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